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Abstract 

Infrared surveys provide essential insights on galaxy evolution. If near-IR studies suggest mild 
evolution of stellar populations with cosmic time, indications of a substantial evolution have 
been seen in the far-IR, although the available information is largely insufficient to delineate 
precise evolutionary properties. 

A consistent picture encompassing all the currently available data may be obtained assuming 
that dust extinction hides the early evolutionary phases of spheroidal galaxies in the optical 
band, while the corresponding dust re-radiation in the far-IR may have been, during the early 
evolutionary phases, orders of magnitude larger than today. Hyperluminous IRAS galaxies might 
be extreme examples of this situation. Additional indications that at least some spheroidal 
galaxies may have been very dusty during their early evolution are provided by recent data on 
high redshift radio galaxies and quasars. 

Galaxies are the likely dominant contributors to the IR background. However, in the frame- 
work of unified models for Active Galactic Nuclei, a large number of nuclei hidden by a dusty 
torus may be expected. Implications for the IR background are discussed. 
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1 Galaxy evolution in the optical-IR. General observational trends 



Optical surveys have constituted for many years the primary tool to investigate the photometric 
evolution of galaxies, reflecting the evolution of stellar populations. CCD detectors have allowed 
to reach extremely faint magnitudes, down to ~ 27-28, i.e. to surface densities as high as 
~ 3 X 10^ galaxies deg~^. A remarkable excess over a smooth extrapolation from brighter 
magnitudes was observed in the B-band counts, beginning at i? ~ 20 and continuing to the 
faintest levels. 

Counts and colors indicated a substantial increase with look-back time of the surface density 
of blue galaxies, consistent with evolutionary stellar population synthesis models provided that 
the universe has a low density {qo ~ 0.05) and galaxies form at high redshifts (zfor — 30). 

The observed redshift distributions of galaxies down to i3 ~ 22-24, however, cannot be 
accounted for by pure luminosity evolution models, at least if the local luminosity function keeps 
relatively flat at its faint end. Various explanations have been proposed, resorting to either 
density evolution ||l| or to luminosity dependent evolution 0, or to a non-zero cosmological 
constant |^]; on the other hand, Gronwall and Koo obtained good fits to the data with 
a standard photometric evolution model with a SMC extinction law, assuming a strong excess 
(still consistent with the data) of low luminosity galaxies over the Schechter luminosity function. 

Crucial constraints, concerning in particular the effects of dust in high-z galaxies [Q], come 
from observations at infrared wavelengths. A great observational stride occurred thanks to the 
impressively fast advances in infrared array technology, which have already made possible to 
reach a limiting magnitude K = 23 fs], |6|, i.e. surface densities ~ 2 x 10^ deg^^ close to those of 
the deepest optical surveys. 

Unlike optical counts, which are highly responsive to young stellar populations, hence prone 
to very uncertain K-corrections and biased against early-type galaxies, near-infrared counts are 
much less sensitive to the detailed star-formation history. By allowing a more straightforward 
interpretation, they are going to set much tighter constraints on galaxy evolution. 

The extremely interesting new result was that counts in the K band show a remarkable 
flattening below K ~ 20, at variance with the steep counts observed in the optical. 

In the far-IR, the first indication that the IRAS galaxy population evolves significantly came 
from 60 fim counts performed at the faintest flux levels ||^, ^ and at medium sensitivities . The 
more detailed information made available by redshift surveys ||l^ suggests a strongly evolving 
luminosity function, with a luminosity evolution rate comparable to that of quasars: 

L{z) = Lo{l + zf^. (1) 

The redshift space, however, was sampled by the IRAS survey only to z ~ 0.2, not enough to 
distinguish between luminosity and density evolution. The limited sampled volume also implies 
that the estimate of the evolution rate in eq. (1) is possibly affected by large-scale structure. 
An underdensity by a factor of roughly two in the local universe has been invoked to interpret 
the very steep counts of galaxies in the optical at B < 18, in particular from the APM survey 



in the southern sky [11 1. 



Near- and far-IR selections emphasize quite different physical processes in galaxies: photo- 
spheric emission of red giant stars, mostly unaffected by diffuse dust, in the K band, and dust 
re-radiation powered by young stellar populations at far-IR wavelengths. To best exploit the 
information content of the two bands, we need a tool to relate the properties of the ISM (gas 
fraction, metal abundance and dust content) with those of the stellar populations. The next 
chapter is devoted to discuss such a tool. 
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2 A comprehensive view on galaxy evolution: the basic recipes 



Evolutionary population synthesis (EPS), models which account for dust effects, have been 
worked out by |l^. Given their large spectral coverage, from ultraviolet to 1 mm, these 
models allow a synoptic view of data relevant to understand galaxy evolution. They deal, in 
a self-consistent way, with chemical and photometric evolution of a galaxy over four decades 
in wavelength. The synthetic spectral energy distribution (SED) incorporates stellar emission, 
internal extinction and re-emission by dust. The stellar contribution, including all evolutionary 
phases of stars born with different metallicities, extends up to 25 fim. Dust includes a cold 
component, heated by the general radiation field, and a warm component associated with HII 
regions. Emission from policyclic aromatic hydrocarbon molecules (PAH) and from circumstellar 
dust shells are also taken into account. 

The star formation rate (SFR) and the initial mass function (IMF) are the basic input 
functions. The most critical parameter, whose variation is able to account for the whole Hubble 
sequence, is the initial value of the SFR, ipQ. Low values of tpQ correspond to a SFR slowly 
decreasing with time, the standard scenario for late-type systems; high values of ■00 provide a 
drastically different behaviour, due to the fast decreasing of the gas fraction, hence of the SFR: 
this corresponds to early-type systems. 

Models successfully match the SED of local galaxies of different morphological type over 
four decades in frequency ([12, |l^). In the following we will refer to models computed with 



a Salpeter IMF and a SFR depending on a power law of the mass of gas: SFR cx M^^^, with 
n = 0.5. This corresponds to a quick metal enrichment and a prompt contamination by dust 
of the interstellar medium: early-type systems become bright far-FIR sources given the high 
optical depth provided by these models during their initial evolution. 

3 Observational tests of the evolution 

The evolution with cosmic time of stellar populations in galaxies can be directly investigated by 
looking at the spectral energy distribution (SED) of galaxies at different redshifts. "Normal" 
high-z galaxies, however, proved to be remarkably elusive: none of the 100 (out of a total of 
104) galaxies with measured redshift in the sample of Colless et al. complete to B = 22.5, 
has z > 0.7; the highest redshift in the complete S < 24 sample of Cowie et al. [15| is 2 = 0.73; 



the measured redshifts of near-IR selected galaxies down to = 20 (nearly 100% complete 
at X < 18 and ~ 70% complete at K = 19-20) are ^ 1, with the exception of one object at 
z = 2.35. 

So far, the only effective method to find high-z galaxies has been optical identifications of ra- 
dio sources [^; far-IR/sub-mm observations of these objects are providing extremely important 
indications on the properties of young galaxies. 

In the following, we will briefly review recent observational results in the near-IR and far- 
IR/sub-mm bands. 

3.1 Near-IR-selected galaxy samples 

Impressive results have been recently obtained by the If A team at Hawaii |16|, [18|. Figures 
1 and 2 compare the data of [16| with synthetic colours yielded by models described in the 
previous section. 

We note that evolutionary effects in the spectra of early-type systems cannot be neglected 
at z > 1; the use of template SED's of local objects to infer photometric redshifts of galaxies 
lacking spectroscopy may be misleading. We have compared, in particular, the colours of the 



16 unidentified galaxies in the K-band selected sample by [16| and found that most of them are 
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likely to be at z > 1. A higher SFR implies bluer colours at z > 1 than expected by simply 
reshifting local SED's. 

Note, in any case, that both the red peak in the colours at z ~ 1 and the subsequent blueing, 
strongly depend on the adopted cosmological model; in particular, colours (I — K) > 5 are 
difficult to obtain in a closed world model. 

Indications in favour of a merging-driven evolution, in which galaxies become more numerous 



but less luminous in the past, have been found e.g. by [18|. Alternatively, we find that, taking 
into account the combined spectral effects of the evolving stellar populations and of extinction 
by dust, all the basic statistical properties of K-band selected galaxies can be explained by pure 
luminosity evolution of the observed local luminosity functions. We account, in particular, for 
the blueing trend observed at > 19 and z ~ 1 as an evolutionary effect due to the higher 
SFR, rather than a change in the galaxy population mixture. 

A world model with qq = 0.15 and Zform = 7 for early-type and ^^form — 2 for late- type 
galaxies is indicated by fits of the observed counts, local luminosity functions, redshift and 
colour distributions. 

3.2 Far-IR and radio-selected samples 

Early-type models including a dusty phase during their initial evolution also match successfully 
the observed spectrum of the ultraluminous galaxy IRAS F10214 + 4724 at z=2.29 ||l9[| . In 
Fig. 3 we compare the observed SED with model predictions; a non-thermal contribution pO| ] 
accounting for 60% of the rest-frame flux at A = 0.1 /im, has also been included. 

Fig. 4 (panels a-c) compares, in the rest-frame, the data on three distant radio galaxies 



having recent sub-mm measurements or upper limits [21, |2^ with our synthetic SEDs. The 
available data can be fully accounted for by "opaque" models like those used by [^] to fit the 
spectrum of IRAS F10214 -|- 4724. The inferred galactic ages are < 0.1 Gyr, thus alleviating 
the problem of galaxy ages uncomfortably high in comparison with the age of the universe as 
well as the difficulty to understand the alignement between optical/infrared continua and the 
(presumably short lived) radio jets |23]. 

4 IR emission of AGN's 

The infrared luminosity is often a significant fraction of the bolometric luminosity of AGN's. 
In the case of radio quiet AGN's, the origin of the IR emission is still debated, although in 
the last few years observational evidence has accumulated in favour of dust reprocessing of the 
primary optical-UV emission. Major issues in favour of dust emission are (cf. and references 
therein): (i) the observed SED's showing a steep rise of the submillimetric continuum between 
1000 and 100 //m which is far too steep for any current non-thermal model, and a local minimum 
in vFi, at A ~ l//m, which is naturally explained by models wherein dust radiates with a limiting 
sublimation temperature Tg ~ 1500 K; (ii) the IR variability and its relationship with UV and 
optical variability with delays among different bands which are easily interpreted in dust models; 
and (iii) the observed low level of polarization. 

Presence of optically thick dust with azimuthal symmetry has been also advocated to ex- 
plain the observed difference between broad- and narrow-line AGN in the framework of unified 



schemes (see |24] for a comprehensive review). As is well known, spectropolarimetric data have 
shown that broad lines also are visible in polarized light in narrow-line nuclei. More recently, 
infrared spectroscopy in the near-IR revealed the presence of broad components in hydrogen 
IR lines (see e.g. ||2^), suggestive of the possibility of a direct view of broad-line regions in 
the IR. Less direct evidence of dust around active nuclei has been found through studies of 
ionization cones in several Seyfert 2 galaxies, which demonstrate that the gas in the host galax- 
ies sees anisotropic ionizing sources. Also X-ray observations revealed high absorbing columns 
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{Nhi ~ 10^^- lO^^cm^^) in many Seyfert 2 galaxies, whereas Seyfert 1 X-ray spectra exhibit 
moderate to low absorption. 

Structure and emission of dust tori or warped discs around AGN's have been investigated 
by several authors (p6|, pO|). In the case of tori, dust is assumed to have standard galactic 
composition, size distribution and sublimation temperature. Various structures have been in- 
vestigated depending on equatorial optical depth, covering factor, shape, radial and axial dust 
density distribution. The models can be broadly divided in two classes. A first one includes 
models with extremely high optical depths {Ay > 500) and extremely compact structure (the 
ratio of the outer rout to the inner radius of the dust distribution ranges from few to several 
tens). Typical problems of these models are that torus emission does not fit the observed spectra 
and surface brightnesses of Seyfert 1 and 2 unless additional dust components are invoked (see 
e.g. 11). 

The second family of models requires less extreme optical depth (Ay ^ 100) and more 
extended tori (nn/rout ~ few hundreds) |Q. These models have the nice properties of fitting 
the observed SED's and of predicting a size of the region emitting at mid-IR extending over 
several tens of parsecs in agreement with the observations. Moreover the predicted anisotropy 
of the 10;um emission is within the limits required by the comparison of the luminosity functions 
of type 1 and type 2 Seyferts 



If the latter models hold, the 12 /im luminosity functions of the two Seyfert types |28| and the 
statistics of the ratio of hard X-ray to 12 /im fiuxes |2^] indicate that not many hidden AGN's 
exist. At most 20% of Seyfert 1 and 50% of Seyfert 2 galaxies could have been misidentified 
as non-active. Moreover the non-linear relation L12 oc L^'j^ suggests that the covering factor 
decreases with increasing source power |2S]. 



5 Contributions to the background flux 

In the absence of dust extinction, the diffuse radiation generated by early stellar nucleosynthesis 
is expected to peak in the near-IR and to have a bolometric intensity: 

O.OOTcV/^ _2 -1 -1 

4ol = ./-,, r erg cm s sr , (2) 

47r(l + zp) 

where ^/ is the galaxy formation redshift, p is the density of the processed material to produce 
a mass fraction of metals. Dust eventually re-radiates this energy at longer wavelengths. 

Estimates of the IR background intensity have been attempted in the near-IR and sub-mm 
bands, where the bright foreground emissions from the Galaxy and Inter-Planetary Dust (IPD) 
are at a minimum. In any case, careful modelling and subtraction of foregrounds is necessary. 
A summary of recent results is shown in Fig. 5. 



5.1 The near-IR background 

The thick solid line labelled "d" in Fig. 5 shows the contribution of starlight as implied by 
models described in §3.1. This estimate is quite a robust one, constrained as it is by existing 
deep counts of galaxies in K band. The predicted spectral shape is also well constrained up to 
A ~ 5;um by our present understanding of photometric evolution. 

The tightest observational limits come from observations of very high-energy (TeV) 7-rays 
in the spectra of the BL Lac object Mkn 421 and follow from the fact that a fiux of extragalactic 
7-rays is attenuated by interactions with ambient photons leading to the production of electron- 
positron pairs (pO[, |31]). These limits are consistent with the near-IR background mostly 



originating from ordinary stellar processes in galaxies. 

Also reported in Figure 5 is the estimated final sensitivity of DIRBE on COBE, which, 
in principle, would allow accurate determinations of the extragalactic background flux from 1 
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to 200 fim. As already mentioned, however, the main hmiting factor are uncertainties in the 
estimates of the foreground components. 



5.2 The far-IR/sub-mm background 

In Fig. 5, the hues labelled "a", "b", "c" show predictions of detailed evolutionary models, while 
the thin line is based on simpler kinematical prescriptions fitting the 60 fim statistic [^]. We see 
an appreciable latitude in such predictions. Note, however, that the extreme curves "a" and "c" 
bound a sort of conservative "allowed region", the former curve, assuming no evolution, being 
already inconsistent with 60 ^m counts, the latter corresponding to an evolution rate probably 
inconsistent with the observed lack of high-z galaxies among faint IRAS sources]^. 

The prediction of the kinematic model implies an higher flux level at mm wavelenghts as it 
assumes a dust temperature distribution constant with time. Physical models predict a lower 
mm background because dust gets warmer in the past, following the increase of the SFR and of 
the average radiation field in galaxies. 

The various observational limits, set by re-analyses of IRAS survey data [10|, by the COBE- 



FIRAS all-sky spectra deconvolved from the Galaxy emission [33] and by rocket observations 



[ 34 1 , are within a factor of a few of the expected extragalactic fiux. The latter turns out to be 



very close to the sensitivity limits of current observations. 
5.3 A mid-IR background from hot dust in AGNs? 

At the minimum between the predicted starlight background in the near-IR and the dust re- 
radiation feature in the far-IR, there is a spectral region (5 < A < 20 /um) in which properties 
of even local galaxies are very poorly known (essentially because of the poor sensitivity of 
IRAS at these wavelengths). Existing data suggest, in any case, that the integrated emission of 
normal galaxies should keep below 3 x 10~^ ergcm"^ s~^ sr~^, unless dust with implausibly high 
temperature characterizes the early evolution phases. 

As previously discussed, dust hot enough to emit at these wavelengths is likely to be present 
in AGNs. The corresponding contribution to the mid-IR background is set by the convolution 
of the local IR volume emissivity of AGNs with their evolution rate. We report in Fig. 5 (line 
marked AGN) a crude estimate based on the 12 //m luminosity function of type 1 and 2 Seyfert 
galaxies by |^] and an evolution rate equal to that of optical QSOs (eq.[l]), extrapolated to 
-^max = 4 with qo = 0.05. The adopted IR spectra for the two AGN types are the average SED 



of Seyfert 1 nuclei and the SED of NGC 1068 |2g. The IR luminosity function of Seyferts 2 has 



been renormalized upwards by a factor 4 to account for the fraction of optically undetected, dust 
extinguished AGNs predicted by the unified AGN model. We see that, with these somewhat 
extreme ingredients, the contribution of dusty AGNs in the range 10-50 fim may exceed those 
of other galaxy populations. 
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Figure captions 



Fig. 1 Colour-redshift relationships predicted by piiotometric evolution models described in 
§2 for galaxies of different morphological types assumed to start their star formation at 
the redshifts Zf indicated, compared with data of |16|. The open symbols show the colours 
of galaxies lacking redshift measurements, plotted at a nominal z = 5. 

Fig. 2 Colour-colour plots for intermediate redshift (0.5 < ^ < 1) galaxies (filled circles) and 
of candidate high-z galaxies (open circles) in the sample of . The lines show predictions 
of the same models as in Fig. 1. The circle and the star mark the synthetic colours at 
z = 1 and z = 2, respectively, of galaxies of each morphological type. 

Fig. 3 Spectral energy distribution of the hyperluminous galaxy IRAS F10214 + 4724 at z = 



2.29 (see [19| for details and references). A non-thermal contribution having a spectrum 
equal to the mean spectrum of Seyfert 1 nuclei [^0| and accounting for 60% of the rest- 
frame flux at A = 0.1 /im, has been included. 

Fig. 4 Fits to the SED's of three high-z radiogalaxies. The ages of the models are specified in 
each panels. The lower thin curves show the assumed nuclear contributions. 

Fig. 5 Extragalactic background light at IR to sub-mm wavelengths. Curve a corresponds to 
no evolution, curves h and c to moderate or strong extinction, respectively, of early type 
galaxies during early evolutionary phases; curve d is the integrated starlight of distant 
galaxies in the near-IR, derived from models fitting the deep K-band counts. Curve S 
shows the contribution of stars in our own Galaxy, at high galactic latitudes. The curve 
labelled "AGN" shows the somewhat extreme estimate of the possible AGN contribution 
to the background described in §5.2. The DIRBE sensitivity as function of A is also shown. 
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